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Abstract. A proper treatment of the non-equilibrium dust formation process is crucial in models of AGB star winds. In this 
paper the micro-physics of this process is treated in detail, with an emphasis on the effects of drift (drift models). We summarize 
the description of the dust formation process and make a few additions to previous work. A detailed study shows that different 
growth species dominate the grain growth rates at different drift velocities. The new models show that the net effect of drift is to 
significantly increase the amounts of dust, seemingly without affecting the mean wind properties, such as e.g., the mass loss rate. 
In some cases there is several times more dust in drift models, compared to the values in the corresponding non-drift models. 
We study the formation of a dust shell in the inner parts of the wind and find that drift plays an active role in accumulating dust 
to certain narrow regions. In view of the results presented here it is questionable if drift - under the current assumptions - can 
be ignored in the grain growth rates. 
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1. Introduction 

A crucial aspect in wind models of asymptotic giant branch 
(AGB) stars is the description of dust formation. It is the radia- 
tion pressure on dust that with the support of the enclosed pul- 
sating star is believed to form and drive the most massive winds 
in these stars. A detailed time-dependent treatment of the dust 
formation is necessary since it occurs far from equilibrium. 

In the second article in this series, Sandin & Hofner (2003b, 
henceforth Paper II), we carried out a thorough study of the 
effects of grain drift on the average outflow properties of sev- 
eral types of time-dependent wind models. The results of wind 
models allowing drift (drift models) were compared with the 
respective non-drift models. A main finding was that drift, in 
most cases, modifies the wind structure to a significant degree 
concerning outflow properties and their temporal variability. In 
particular wind models that use a more realistic gas opacity 
are affected. The work presented in this article is based on the 
model description given in Paper II. Grain drift has so far been 
dynamically included through the use of a separate equation of 
motion for the dust. It was, however, not included in the pro- 
cesses describing dust formation. In regard of the drift induced 
changes found in previous results it is questionable if this treat- 
ment is adequate. 
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In this article we caiTy out a closer study of the micro- 
physics of dust formation in the wind forming region of time- 
dependent models, allowing drift. We begin by modifying the 
description of the grain growth process. Thereafter, results are 
discussed to assess the role of drift to the wind formation - and 
the formation of dust shells. The purpose of this article is to 
focus on the understanding of the wind formation using a few 
typical models. A study closely related to this paper by Kriiger 
& Sedlmayr (1997, henceforth KS97) was concerned with the 
grain size distribution and dust formation in stationary mod- 
els including drift; their conclusions are different form those 
found here. The detailed treatment of the dust material prop- 
erties in time-dependent wind models was recently addressed 
in a study carried out by Andersen et al. (2003, henceforth 
AHG03). The results showed significant differences depending 
on the adopted properties of the dust. 



The modifications we carry out to include drift in dust for- 
mation processes are first described in Sect. 2. Then the mod- 
eling procedure is presented together with a discussion on av- 
eraged outflow properties in Sect. 3. The consequences of al- 
lowing drift for the formation of a dust shell, and the details of 
the micro-physics of dust formation are discussed in Sect. 4; 
followed by the conclusions in Sect. 5. 
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2. Including drift in the dust formation description 2.1. The treatment of grain growth neglecting drift 



Like in the earlier articles in this series we distinguish between 
three interacting physical components of the wind. The com- 
ponents are the gas, the dust, and the radiation field. Each of 
these is described by coupled conservation equations that in- 
clude exchange of mass, energy and momentum between all 
three components. A thorough description of the physical sys- 
tem, the gas-dust interaction, and the numerical method was 
given in Paper I. The effects of stellar pulsations and an im- 
proved treatment of the gas opacity were added in Paper II. 
The work presented here is based on this most recent formu- 
lation. New in this article is the improved description of dust 
formation, which is extended to include effects of drift. 

The dust component is assumed to consist of spherical par- 
ticles made of amorphous carbon. Dust formation is described 
using the so-called moment method (Gail et al. 1984, Gail & 
Sedlmayr 1988, and Gauger et al. 1990 henceforth GGS90), 
involving four moment equations. The moments represent cer- 
tain (average) properties of the grain size distribution function 
(see e.g. Paper II, Table I). Assumptions of this method are that 
grains are large enough that their thermodynamic properties do 
not depend on the grain size, and only molecules with a few 
monomers contribute significantly to the growth process. The 
moment equations can be written as (cf. GGS90), 



-Kq + V-{K^v) = J 
ot 



d_ 



n 1 



-K, + V • (K„v) = - -K„.i + N^'''j : (l<n<d). (2) 



Here d is 3 for spherical grains, Ni is the lower size-limit of 
macroscopic grains (we use M = 1000 carbon atoms), and 
V is the dust velocity. Equation (I) describes the grain pro- 
duction, and the other three equations different properties of 
grains. In the form these equations are presented ^ is the grain 
nucleation rate. Nucleation describes the interchange of par- 
ticles between the dust and gas phases (i.e. the formation of 
seed particles). Like in previous work we adopt a stationary 
description (cf. Gail & Sedlmayr 1988, Sect. 4), using clas- 
sical nucleation theory. We are, however, aware that this is a 
problematic point and that this assumption introduces quantita- 
tive uncertainties in the models. For a recent discussion on the 
nucleation rate and issues related to currently available descrip- 
tions see, e.g., AHG03 and references therein. I/t denotes the 
net grain growth rate (see next subsection). 

Like in several of the earlier wind model articles we adopt 
a C-rich equilibrium chemistry. The atomic and molecular 
species involved in the dust formation are H, H2, C, C2, C2H, 
and C2H2, where the last four contribute to grain formation pro- 
cesses. Dust formation has so far been treated accounting for 
nucleation, homogeneous growth, thermal evaporation, chem- 
ical growth, and chemical sputtering (see the following sub- 
section). One assumption used in the description of these pro- 
cesses is that drift between gas and dust is negligible. In this 
study drift is the key feature and in the following subsections 
we discuss its implementation. In particular the net growth rate, 
T"\ is affected. 



In this subsection we describe grain growth without drift in or- 
der to identify the terms that need to be modified when allowing 
drift. GGS90 wrote the grain growth rate in the following form, 



1 



1 



1 



(3) 
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PC ^gr,h ^ev,h ^gr,c ^sp,c 

The first term, 1 /Tgr,h, describes homogeneous growth by addi- 
tion of carbon i-mers (C, ; monomers). The second term, 1 /Tev,h, 
represents thermal evaporation from the surface of grains; as 
such it only depends on the properties of the dust. Reactions 
involving the two molecules C2H and C2H2 are described by 
chemical growth (also called heterogeneous growth) 1/Tgr,c and 
chemical sputtering 1 / Tsp,c . The physical form of these terms is, 
when chemical equiUbrium and a negUgible drift are assumed, 
as follows. 
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(7) 



The summation over i labels the / different carbon j-mers ac- 
counted for, in this case C and C2. Likewise m labels the M, 
different reactions for molecules involving /' carbon atoms, in 
this case C2H and C2H2 for / = 2. Moreover the gas and the 
dust temperatures are specified by Tg and r^. TC,, 'Ki m, and "TC,"^ 
are the dissociation constants that can be used to calculate par- 
tial pressures of the relevant molecules if chemical equilibrium 
is assumed (cf. GGS90 and references therein). The average 
sticking and reaction efficiencies are given by a, and a)^^, re- 
spectively (see Sect. 2.2.1). S denotes the super- saturation ra- 
tio, 



S 



(8) 



^C,sat(rd) 

where PciTg) is the actual partial pressure of carbon atoms in 
the gas phase. The denominator, Pc,sit(Td), is the vapor (satu- 
ration) pressure of carbon atoms over a solid carbon surface, 
specified at the dust temperature T^. Like in earUer work we 
use (Gail & Sedlmayr 1988, Sect. 7), 

86300 

logPc,sat = ^ + 32.89 (9) 

which is derived for graphite and is valid in the temperature 
regime 500K < T < 1500K. The (hypothetical) monomer 
surface area Aj is given by. 



Al = 4nr, 



,2^4^(3:?lcm„ 



\ 4«-pgr 



(10) 
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where ro is the monomer radius, JIq is the atomic weight of a 
carbon atom, m„ is the atomic mass constant, and pgr the mass 
density of the condensed grain material. As in the previous ar- 
ticles we use the value for graphite, pgr = 2.25 gcm"^. Note, 
however, that a value for amorphous carbon would be more 
consistent (cf. AHG03). The number density of /-mers and the 
individual molecules involved in the building of dust grains are 
denoted by /(«') and n^m, respectively. The average thermal ve- 
locity of carbon /-mers, radicals, and molecules in the gas phase 
is given by the usual Maxwell-Boltzmann mean Vi(,m)' for each 
species. Particles that hit a grain surface cannot have a velocity 
directed such that they are moving away from the dust grain; 
it is easily seen that the integrated infalling flux of particles 
through a (planar) surface is given by O = /(Ovi/4. For our 
purposes the velocity in Eqs. (4)-(7) is defined as. 



Vi(,m) 



1 / SfeB^g 

4 \nJli(,„)ma 



(11) 



where yi,(,,„) is the atomic weight of the particle species in ques- 
tion (see Table 1). 

2.2. The treatment of grain growth including drift 

The growth rates are modified in several ways when a relative 
motion between the gas and dust phases is allowed. Before we 
sum up the new rates in Sect. 2.2.3 we discuss the necessary 
physical difi'erences. Related works on stationary winds were 
carried out by Dominik et al. (1989) and KS97 (Sect. 2, who 
included fewer processes in the grain growth rates). 

2.2.1. The drift-dependent relative velocity and 
sticking coefficients 

The approximate relative velocity between particles moving 
with both a thermal velocity and a drifting velocity can be 
written as a root-mean- square of the sum of these two veloc- 
ities (Draine 1980, Sect. 5a). Since dust grains are much larger 
compared to gas particles, we again define our relative velocity 
by the same argument leading to Eq. (11), 



1 



>Vi(,m) 



vpia) 
16 



(12) 



where VD(fl) = v(a) - m is the grain size dependent drift velocity, 
with a denoting the grain radius. A drifting motion increases 
the flux of particles hitting a grain and thereby enhances the 
formation processes (see below). 

Since we treat the drift velocity dependence of sticking co- 
efficients (or,) and of reaction efficiencies (a'r ) in the same 
way, we in the following refer to both these quantities as stick- 
ing coefficients, a,(,m)- Gaseous (hydro)carbon species may 
bond to radical surface sites on dust grains (chemisorption). 
The binding energy of the (hydro)carbon to the surface, E},, de- 
pends on the (unknown) surface morphology of the particle in 
question. The translational energy of the hitting particle must 



Fig. 1. The sticking coefficient a (Eq.(13)) drawn as a function of 
the drift velocity for each of the four hydrocarbon species (up- 
per panel). These species that are all involved in the grain growth 
rate chemistry are: C (dash-dot-dotted line), C2 (dash-dotted line), 
C2H (dotted line), and C2H2 (solid line). In addition the new/old ra- 
tio of the term wa is shown in the lower panel. The figure shows that 
a(C2H2) drops at a much lower drift velocity compared to the other 
three species, also note the steeply increasing grain formation effi- 
ciency (for all species) even at low drift velocities. The relative ve- 
locity w is calculated using a gas temperature of Tg = 1200 K. Note, 
however, that w is mostly dependent on the drift velocity when this is 
high. 



be adsorbed by the target, or it will bounce off. As a veloc- 
ity dependent expression for the sticking coefficients we use 
the relation given by KS97 (note that the exponential form of 
this relation is based on a study involving the sticking of argon 
atoms onto argon-covered ruthenium). 



ai(,m) = ac,i(,m) exp 



/'0.5j?l,(,m)muw|_, 



,m) 



4£b,,-, 



K,m) 



3\ 



(13) 



' The term x^^m) is a short notation referring to both x, and x, 



The nominator of the exponent in Eq. (13) specifies the relative 

kinetic energy of the particle hitting a dust grain. In this expres- 
sion it is assumed that a translational energy 4 times the binding 
energy can be adsorbed by a dust grain. We use binding ener- 
gies based on quantum-mechanical calculations of the associa- 
tion of carbon species to sp^ diamond surfaces. Binding ener- 
gies and particle masses for each species are given in Table 1. 

In contrast to the approach chosen by KS97 we adopt con- 
stant sticking coefficients ac,i(,m) (smaller than one) in front 
of the exponential term in Eq. (13). This is done to keep the 
same sticking coefficients as in previous models in the limit of 
zero drift velocities. For comparison reasons we adopt the same 
numbers we have used in earlier articles (e.g.. Paper I, Paper II, 
and Hofner et al. 1995 to mention a few), cf. Gail et al. (1984) 
and Table 1 (these are the same sticking and reaction efficien- 
cies used in Eqs. (4)-(7)). Kriiger et al. (1996) and KS97 study 
sticking coefficients in more detail. The sticking coefficients as 
given by Eq. (13) are plotted as a function of the drift velocity 
for each used hydrocarbon species in Fig. 1 . Note that the stick- 
ing probability of C2H2 drops much faster with increasing drift 
velocity than those of the other species do; it quickly drops to 
zero when vr,{a) exceeds about lOkms"^ The sticking proba- 
bility of the other three species are more or less unaffected even 
when VD(a) = 40 km s"' . 

The flux of gas particles hitting a grain surface increases 
with the drift velocity (Eq. (12)), thereby increasing the growth 
rate. If the drift velocity is too large, gas particles are too ener- 
getic to stick to a grain surface and instead bounce off, inhibit- 
ing further growth. However, the situation is ambiguous since 
different species contributing to the grain growth have differ- 
ent sticking probabilities. The increasing efficiency of grain 
growth/destruction with the drift velocity is illustrated using 
flie term H';(,m)0'i(,m)/v,(,m)Q'c,!(,m) in Fig. 1 (lower panel). Note 
that the grain growth efficiency increases fast with the drift ve- 
locity; it is about two times larger already for VD(a) ~ 2kms'"' 
(more dust is indeed formed in the new drift models. Sect. 3.3). 
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Table 1. Molecular properties of the (hydro)carbon species taking part 
in the grain growth process. From the left the columns give: the m 
used in the growth rates; the molecule species; the atomic weight; the 
(constant) sticking coefficient (Eq. (13); these are the same numbers 
used for a; and in Eqs. (4-7)); the binding energy for each species 
to a sp^ carbon surface; and the corresponding source of binding en- 
ergies in the literature. 



i,m 


mol. 


A 






references for E^, 






[Wu] 




[eV] 




1 


C 


12.011 


0.37 


4.0 


priv. comm. . . . 


2 


C2 


24.022 


0.34 


5.0 


... K. Larsson (2003) 


2,1 


C2H 


25.030 


0.34 


6.5 


Larsson et al. (1993) 


2,2 


C2H2 


26.038 


0.34 


0.27 


Larsson et al. (1993) 



2.2.2. Non-thermal sputtering; the erosion of dust 
grains by energetic abundant gas particles 

A process unique for situations involving drift is non-thermal 
sputtering in which gas particles (projectiles) at high drift ve- 
locities are energetic enough to tear off carbon atoms from a 
grain surface (target). The corresponding rate of grain destruc- 
tion is given by KS97 as, 



— = nrlvD(a)y]njYspjiEj) 

Tsp,n . 



(14) 



where the sum over j is taken for the different contributing par- 
ticle species in the gas. The sputtering yield for each species is 
denoted by ^spj- The number density and the kinetic energy of 
the projectile gas particles are given by nj and Ej, respectively. 
Ej is simply, 



(15) 



where Jlj is the atomic weight of the projectile particle. In 
the case of a solar chemical composition of the gas, contri- 
butions from other elements than hydrogen and helium can 
be neglected (Woitke et al. 1993, see Fig. 1). Furthermore, 
non-thermal sputtering is more or less neghgible at drift veloc- 
ities below 30 km s"' . The drift velocity has earlier (Paper I,II) 
been found mostly to be lower than this. Nevertheless, non- 
thermal sputtering is included here since we want to get an 
idea of its relevance, if any, to the wind structure. For the 
yields in Eq. (14) we adopt the empirical expressions given 
by Bohdansky et al. (1980, 1981, also see Woitke et al. 1993 
and the additional references given therein). 



f 0.0064. 

i -.EjK E^j 



Here jj is given by. 



(J?lc + Mj) 



2 ■ 



(16) 



(17) 



Moreover is the atomic weight of the projectile particles 
(i.e. about 1 and 4 for H and He, respectively). £thj is the sput- 
tering energy threshold and determines when sputtering starts 



to occur, it is specified as, 

^b,C 



r/1 - 7j) 



^>0.3 



(18) 



Eb,c(= 4.0 eV) denotes the surface binding energy of individual 
carbon atoms on a dust grain. An alternative to Eq. (16) is given 
in the semi-analytical yields presented by Tielens et al. (1994). 

2.2.3. The growth rates accounting for drift 

All growth rates given in Eqs. (4)-(7) describe the interaction 
of gas particles with dust grains. At the presence of a non-zero 
(average) drift velocity between these particles sticking coef- 
ficients and fluxes of particles hitting grain surfaces are modi- 
fied. Thereby the homogeneous and chemical growth rates and 
the chemical sputtering are affected. The thermal evaporation 
(Eq. (5)) is, however, not affected as it only depends on the 
properties of the dust. Allowing drift we replace the total grain 
formation rate Eq. (3) with. 



1^^^ L + J 1 

^ ^gr,h T"ev,h ^gr,c ^sp,c ^! 



1 

sp,n 



1 



1 



' sp,n 



(19) 



where Tq' is the total drift-dependent growth rate excluding 
non-thermal sputtering. The growth rates given in Eqs. (4)-(7) 
are replaced with. 



1 

1 

^ev,h 

1 



'ViAiaiiwd-WifiUt) 



V -A - fr A 1 ^(^d) 
> lAxOi ■ Vifii, t) ■ — — 



1=1 



M, 



(20) 



(21) 



(22) 



1=1 



m=l 
Af, 



l^WTCUTd) 



(23) 



A fundamental difference, when compared to the earlier 
formulation, is that the rates now are grain size dependent 
through the drift velocity. When used with the moment equa- 
tions (Eqs. (1) and (2)) this poses a problem since these do 
not readily allow for the use of grain size dependent rates. 
Moreover, our system of equations currently only includes one 
equation of motion for the dust, and thereby one mean grain 
velocity. A binned grain size distribution would require one ad- 
ditional equation of motion for each bin. 

In Paper I we argued that the dust velocity can be de- 
scribed using one mean quantity. With the current physical and 
methodological htnitations we apply this assumption to grain 
growth as well, hence. 



vv>{a) = VD«rd» = vd 



(24) 
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where (r^) is the average grain radius. Our study is hereby lim- 
ited to qualitative effects of drift-dependent grain growth, and 
we do not study detailed effects associated with the grain size 
distribution. 

3. Modeling procedure and general results 

3.1. Modeling procedure 

The system of equations is discretized in the volume-integrated 
conservation form on a staggered mesh. The equations are 
solved implicitly using a Newton-Raphson algorithm where the 
Jacobian of the system is inverted by the Henyey method (cf. 
Paper I). In difference to the treatment of models calculated 
earlier we adopt a length scale of the artificial viscosity and 
artificial mass diffusion twice as large compared to the values 
used so far (/ = = 7.0 X 10"^ see Eqs. (13) & (14) in 
Paper I). In some cases this higher value, in regions of very 
low dust densities, reduces large unwanted radial variations in 
the drift velocity (cf. Paper I, Sects. 3.2 & 4.2); apart from that 
it has been found to only marginally affect the wind structure 
as a whole. 

The modeling procedure is as follows. The wind model 
is started from a hydrostatic dust-free initial model where the 
outer boundary is located at about 2/?,. All dust equations are 
switched on at the same time. Dust starts to form whereby an 
outward motion of the dust and the gas is initiated. The ex- 
pansion is followed by the grid to about 25 /?„ where the outer 
boundary is fixed allowing outflow. The drift models evolve 
for about 50-200 P. The low density in the Planck mean mod- 
els results in a 'low mass' envelope that quickly is depleted 
of material. Therefore the time period for which the average 
outflow properties are calculated is shorter than the totally cal- 
culated time of the respective model (cf. Sect. 4.2, Paper II). 
The instants studied in Sect. 4 are selected at times before a 
significant fraction of the envelope is lost. 

To describe the effects of stellar pulsations on the atmo- 
sphere we use a sinusoidal, radially varying inner bound- 
ary, located at about 0.91 (above the region where the k- 
mechanism supposedly originates). An inflow of mass through 
the iimer boundary is not permitted. 

3.2. Selection of model parameters 

As was found in Paper II, wind models show a wide range of 
values in outflow properties and variability. And depending on 
what physical description is used for the gas opacity, the re- 
sulting density structure produces completely different winds. 
With this knowledge it is hard to choose a set of model param- 
eters for one model that will be representative of a much larger 
set of wind models. 

Models calculated with a Planck mean gas opacity result in 
much more realistic density structures compared to the mod- 
els calculated with a constant gas opacity (see e.g. Sect. 2.2 in 
Paper II). The former treatment is adopted here. Moreover, in 
Paper II we found some models to be multi-periodic, meaning 
that the properties at the outer boundary vary with a period that 
is an integer multiple of the stellar piston period. Such a be- 



Table 2. Model parameters, cf. Sect. 3.2. The model names are given 
in Col. 1. Starting from the left the remaining columns specify: the 
stellar luminosity L,; the effective temperature T^g; the pulsation pe- 
riod P; the pulsation amplitude Awp; and the carbon/oxygen ratio 
sc/eo. The stellar mass M, is set to 1.0 Mq in all models. 



model 


L, 




P 


A Mp 


Sc/eo 






[K] 


[d] 


[kms-i] 




P13C16U6 


1.3 X 10"* 


2700 


650 


6 


1.6 


P13C14U6 


1.3 X lO"* 


2700 


650 


6 


1.4 


P10C18U4 


1.0 X lO-* 


2790 


525 


4 


1.8 



havior is not seen in the typical model, and intentionally three 
sets of model parameters are selected that previously resulted 
in irregular wind variability; see Table 2. 

The winds of the two drift models P10C18U4 and 
P13C14U6 in Paper II were both found to give a smaller mass 
loss rate than that of P13C16U6. The latter model can be de- 
scribed as an 'average' model in terms of other properties (see 
Tables 2 & 5 and Sect. 4.2 in Paper II for further details). In the 
following model P13C16U6 is discussed in detail, and the two 
other models are presented for comparison with this model in 
Sect. 3.3. 

The modified physical descriptions of the grain growth 
rates discussed in Sect. 2 each affect the wind structure to dif- 
ferent degrees. We want to assess the importance of each pro- 
cess by itself and therefore discuss four different versions of the 
model. One without drift - a so-called position coupled (PC) 
model - and three drift models. Of the latter three one is cal- 
culated without the modifications in the growth rates (i.e. using 
Eq. (3)). The second and third models adopt the drift modified 
growth rates, and the third in addition includes a description of 
non-thermal sputtering. 

3.3. Differences in averaged outflow properties 
caused by drift in dust formation 

The main issue in Paper II was to study temporal variations 
of outflow properties which result from drift. To compare with 
those findings we here cormnent on outflow properties of three 
different kinds of drift models. The outflow properties (q), with 
the corresponding standard deviations (as), and relative fluctu- 
ation amplitudes (r - os/q), are given in Table 3. The values of 
the models with the suffixes -p— and -<i__ are, with one excep- 
tion, all identical to the corresponding PC model values given 
in Paper II, Table 5. The exception is model P13C16U6-d_ 
which has been recalculated here^. 

The effects of drift-dependent dust formation excluding 
non-thermal sputtering are studied using model P13C16U6- 
dv-. The only apparent difference in the outcome of 
P13C16U6-c/v_ and P13C16U6-t/Ks (which includes non- 
thermal sputtering) is the higher average drift velocity in the 
latter model. This higher value could be a result of how the 

^ The difference is due to the slightly different time-interval used in 
the calculation of the average properties, and the different length scale 
adopted for the artificial viscosity/diffusion. 
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Table 3. Quantities temporally averaged at the outer boundary for models P13C16U6, P13C14U4, and P10C18U4 (see Sect. 3). Physical 
features of the different versions of the model are indicated with the model name in the first column. Drift models are denoted by a 'rf' and PC 
models by a 'p' . Moreover a 'v' indicates the use of the "relative" velocity given in Eq. (12), and an the use of non-thermal sputtering. An 
unused feature is indicated by an The following columns give the mean mass loss rate {M), the mean terminal velocity {u„), the mean degree 
of condensation (famd), the mean dust/gas density ratio {pn/p}, and the mean drift velocity (I'o), respectively. In addition the standard deviation 
(cTs) and the relative fluctuation amplitude r{=crslq) are specified for each quantity (q). The values shown in bold face of the quantities of the 
-dv- and -dvs models indicate that they differ significantly (by > 10%) from the corresponding values of the -d.. models. All models show an 
irregular temporal variability. 
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7.3 (0.33) 
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0.67 
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(1.3) 


0.068 


33 


(32) 
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26 


(93) 


3.5 
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(1.4) 


0.074 


32 
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3.2 
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(2.7) 
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(24) 


1.0 
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(30) 


3.6 


10 
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(1.5) 


0.79 
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(3.4) 


0.40 


39 


(25) 


0.64 


11 


(19) 
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8.4 


P10C18U4-p_. 


1.1 


(0.49) 


0.45 


14 


(0.77) 


0.055 


16 


(2.6) 


0.16 


7.2 


(1.2) 


0.17 




P10C18U4-rf_. 


0.82 


(0.12) 


0.15 


11 


(0.24) 


0.022 


12 


(5.3) 


0.44 


5.4 


(2.9) 


0.54 


4.5 


P10C18U4-rfvi 


0.81 


(0.32) 


0.40 


19 


(0.50) 


0.026 


46 


(21) 


0.46 


21 


(14) 


0.67 


7.0 



mean of the drift velocity is calculated. The higher value possi- 
bly arises due to the higher drift velocities found in the regions 
in front of shocks in model P13C16U6-(iv.s (see Sect. 4.2). 

The values of the new sets of models, i.e. those with a 
suffix -dv- and -dvs, differ significantly (defined as > 10%) 
in most quantities and fluctuation amplitudes from the corre- 
sponding values of the two old sets of models -d— and -p 

In particular, the higher values of both the degree of conden- 
sation and the dust/gas density ratio in all new winds indi- 
cate a more efficient dust formation when drift is accounted 
for, even at low drift velocities of only a few km s~^. That the 
wind acceleration works more efficiently - as a consequence 
of increased amounts of dust (seen in the dust/gas density ra- 
tio) - is indicated in the higher terminal velocity of models 
P10C18U4-^/vw, P13C16U6-t/v_ and Pl3Cl6\J6-dvs. The situ- 
ation appears less certain in P13C14U6-dvs where the value on 
the terminal velocity falls between the values of P13C14U6-p_ 

and P13C14U6-(i As discussed in Paper II, PC model winds 

with a terminal velocity < 10 kms"' have no counterpart in 
drift models - no wind is produced. Model P13C14U6-p_ is 
a border line case, resulting in a much lower mass loss rate 

in P13C14U6-d Due to the lower carbon abundance in this 

model (defined in the carbon/oxygen ratio), there is not enough 
material to form as much dust as in the other -dv.s and -dv_ mod- 
els. In contrast to the case of the terminal velocity the mass loss 
rate appears to be insensitive to the amounts of dust. AH -dvs 
(and -dvS) models show unchanged values when compared to 
the respective -d— models. 

The issue of a larger variability in drift models was dis- 
cussed in Paper II (Sect. 5.2). It was found that the variability in 
drift models mostly is larger than it is in PC models. Similarly 
the relative fluctuation amplitudes in the new -dvs models in 



most quantities differ from the corresponding values of the 
models. Concluding this section we note that the new results 
presented here seem to be difficult to reproduce without actu- 
ally allowing drift in the calculation of dust formation. 



4. Discussion 

In this section we look closer at the region of wind formation. 
First we study the formation of a dust shell in detail in Sect. 4.1. 
Then, in Sect. 4.2, we look closer at the efl'ects on the micro- 
physics in dust formation caused by the addition of drift to the 
grain growth rates. 

4. 1 . The formation of a dust shell allowing drift 

Physical structures of drift model winds generally show a 
larger variability than corresponding PC model winds do (see 
Sect. 3.3, and Paper 11). This behavior is caused by an al- 
lowed dynamic accumulation of dust to narrower regions be- 
hind shocks; here we study how this accumulation might take 
place. The PC model P13C16U6-p_, for instance, shows a 
much smaller variability than the three drift models do and in 
fact reminds of a stationary wind. In the following we study the 
drift model P13C16U6-<ivi, which is calculated including both 
drift-dependent dust formation and non-thermal sputtering. As 
has been pointed out earlier (in Sect. 3) physical structures of 
different wind models show great variations and it is probably 
not safe to generaUze details such as numbers and physical lim- 
its found for this wind model to all other possible cases. This 
study is a complement to the dust shell formation studies ig- 
noring drift given in, e.g., Fleischer et al. (1992, 1995). 
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Fig. 2. Evolutionary sequence of the radial structure of the drift model P13C16U6-(ivi covering five instants of one complete dust formation 
cycle. The first column shows the dust formation at an arbitrarily selected time defined as (fo+)0.00 P (in units of the stellar pulsation period). 
The three following columns are shown at 0.57 P, 1.12 P, and 1.63 P, respectively. The structure of the repeated dust formation cycle at 2.00 P 
is shown by a solid line in the first column, illustrating the aperiodicity of dust formation in the current model. From the top the panels in each 
column show: a) the gas velocity u; b) the gas density p; c) the gas temperature T^; d) the drift velocity Vd; e) the degree of condensation 
/cond; f) the dust/gas density ratio pi/p; and g) the net growth rate r"' . All panels are shown at the same scale in all four columns. The process 
illustrated in this figure shows that a nonzero drift velocity allows for the formation of narrow features in the dust already in the inner parts of 
the wind, cf. Sect. 4.1. 



One arbitrarily selected cycle of the dust shell formation 
is illustrated in Fig. 2. The dust formation is in the following 
discussed in detail for each presented instant. 

1^' column, 0.00 P: One of the shock waves (at 2.4/?,) in- 
duced by the stellar pulsations has reached the region of effi- 
cient grain growth at 1.6 - 2.4/?,, Figs. 2a-b. The carbon in 
ffie gas quickly condenses into and onto grains, forming a new 
dust shell. Fig. 2e; the current maximum degree of condensa- 
tion, /cond ^ 20%, has been reached during 0.23 P (not seen in 
the figure). Small amounts of dust grains are present inwards 
to a second shock at 1.6/?,, inside of which it is too hot for 
grains to exist (and they instead evaporate). Grain growth is 
most efficient in the dense region behind the dust shell, where 
r < 2.4 Rt (Fig. 2g). The present dust heats the region behind 
the dust forming shell, causing a 100 K temperature step (back- 
warming). Fig. 2c. The dust experiences an outwards directed 
radiative pressure originating in the central star, but the dust 
shell is not yet massive enough to push ffie gas outwards by 
itself. The drift velocity is low. Fig. 2d, vd < 3kms"' in all 
of the inner region; the relocation of dust is consequently slow 
as it takes ffie dust about 2.0 P to travel a distance of 1.0/?, 
relative to the gas at this velocity. 

2"'' column, 0.57 f.- The dust shell (now at 3.3/?,) has 
grown and is massive enough to drag ffie gas outwards as a 
consequence of the radiative pressure. A maximum of 80% of 
the available carbon has currently condensed onto grains in ffie 
shell. Another dust shell is forming in ffie region behind ffie 
second shock at 2.0/?,; a region of very efficient grain growth. 
The drift velocity is larger than in the previous column, reach- 
ing about 15kms~^ Note the depletion of dust in the form- 
ing dust shell behind 2.0/?, caused by a drift velocity reaching 
about 5.0 km s"^ 

3 column, 1 . 1 2 P.- The dust in the previously forming dust 
shell (at 2.0/?, in the previous column) has been fully diffused 
into ffie region in front of ffie shock, now at 2.5 /?,. The remains 
are seen as a bump in the degree of condensation at about 4.0 /?, 
(Fig. 2e) behind the dust shell. The physical conditions, in the 
form of a low gas density and a small amount of dust grains, in 
ffie region between 2.5 and 4.8 /?, allow for a large drift veloc- 
ity. Here it reaches values above 30kms"', where non-thermal 
sputtering is active (see Sect. 4.2). Relative to the gas the dust 
moves 2.3/?, in 0.5 P at 30kms"\ and ffiereby allows for a 
quick accumulation to dense regions. This column shows ffie 
same instant that will be discussed in Sect. 4.2. 

4* column, 1.63 P.- The large drift velocity feature has 
moved outwards and almost caught up wiffi ffie dust shell in 
front (now at 6.0/?,). In the process, the dust in front of the 
feature has been "swept up" and is now mostly contained in 



a narrow shell which is merging wiffi ffie original dust shell. 
The drift velocity behind ffie feature, at radii < 4.5/?,, is low. 
A new dust shell is about to form at 2.3/?,, behind a third gas 
shock (emitted 2.00 P after ffie ffist shock discussed in the P' 
column). 

2.00 P (solid line in column 1 ): In this frame two pulsation 
periods have passed since the filled line in ffie first column. The 
radial location of the next forming dust shell, now at 2.8/?,, is 
diff'erent from the one forming at 0.00 P (ffien at 2.4/?,), indi- 
cating a non-periodic dust formation cycle. 

The study of this dust shell formation cycle illustrates how 
dust may be accumulated to ffie regions behind shocks in drift 
models. While the wind is calculated out to 25 P,, our discus- 
sion of the dust formation is limited to the innermost region of 
the wind, inwards of about 7 P,. The outer region is subject to 
complex interactions of boffi physical and numerical character 
that complicates an interpretation. It is unclear if the conclu- 
sions drawn here can be extrapolated to a more extended re- 
gion. 

Next we discuss changes in ffie micro-physics, due to al- 
lowed grain drift in the dust formation, in more detail. 

4.2. The effects of drift on the dust formation process 

The physical conditions in cool C-rich stars are often suitable 
for efficient dust formation in a dynamic region around about 
2.0-2.5 P,. Different atomic and molecular species dominate 
the growth in different temperature (and pressure) regimes. At 
temperaffires Jg < 1500 K it is C2H2, for 1500 K < Tg < 
1900 K C2H, and for Tg > 1900 K free C atoms dominate (see 
Gail & Sedlmayr 1988, Fig. 1). Of these the two former parti- 
cles are ffie dominating carbon bearing species in the wind. 

Figure 3 shows the radial strucffire in ffie inner parts of 
model P13C16U6-Jvs at one instant of ffie dust formation 
cycle. This particular instant is selected as it simultaneously 
shows several distinguished features of the dust formation pro- 
cess (the same instant is shown in the ffiird column in Fig. 2). 

A massive dust shell coinciding wiffi a shock in ffie gas 
has formed and is moving away from ffie star at about 4.8 P, 
(Figs. 3a,c,d; also see the previous subsection). Due to the 
physical conditions, nucleation is currently sharply hmited to 
ffie region between ffie two outermost shocks, where 2.4P, < 
radius < 4.8P,, and it is most efficient in the innermost 
(and most dense) part, Fig. 3g. The same region provides 
suitable conditions for a large drift velocity, which reaches 
30kms"^and above, by a low gas density (Fig. 2b, Col. 3) and 
a small amount of dust grains. Fig. 3b. Few gridpoints are lo- 
cated in the region of a large drift velocity, resulting in a low 
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Fig. 3. Illustration of the effects of drift in the inner parts of the stellar wind of model P13C16U6-(i__ (solid line) and P13C16U6-t?v5 (dotted 
filled line). The panels show: a) the gas velocity u; b) the drift velocity Vd; c) the dust density pi, d) the average grain radius (rd); e) the gas 
temperature T^, f) the chemical growth (without sputtering) Tg/^(C2H) normalized to the net growth rate (excepting non-thermal sputtering) 
Tq'; g) the nucleation rate J,; h) the chemical growth (without sputtering) Tg/j.(C2H2), also normalized (see f); i) the net growth rate excepting 
non-thermal sputtering t^' (Eq. (19)); and j) the non-thermal sputtering t'^^^ (Eq. (14)). All plots are drawn as a function of the stellar radius 
R, (lower axis), alternatively in astronomical units (upper axis). The dots on the contour represent individual grid points, a majority of which 
are located to the shocked regions. The gray horizontal lines are guides. Acetylene is normally the main growth species (h) in the outer cool 
and dilute regions of the wind (modeled here). However, when the drift velocity is larger than about 10kms"'(b) the sticking coefficient a 
for acetylene quickly drops to zero (see Fig. 1), and the radical C2H takes over this role (f). Note that non-thermal sputtering (by helium; j) is 
present when the drift velocity reaches about 30-35 km s"'. 



resolution. The inner boundary of the peak in the drift velocity 
coincides with the shock front in the gas at 2.4/?.. 

The most abundant hydrocarbon molecule in most parts 
of the modeled envelope is C2H2 (acetylene), which also is 
the main growth species in most parts but the region around 
2Rt, where the radical C2H instead dominates, Figs. 3f,h. Very 
small amounts of dust reside in the innermost region (< 1.5/?,, 
Fig. 3c), where the chemistry and resulting properties therefore 
do not play a role. Contributions to the grain growth of the re- 
maining two species, C and C2, are negligible in the current 
model. 

Since the molecule C2H2 compared to, e.g., C2H has a 
lower binding energy to the surface of a dust grain, it is in our 
calculations quickly exchanged as a primary growth element 
by other species when the drift velocity reaches lOkms ' and 
above. In this case by C2H, which sticks to the dust grains even 
if the drift velocity is larger than 40kms"'(see Fig. 1). This is 
nicely illustrated in a comparison of the radial location of the 
peak of the drift velocity in Fig. 3b with the respective growth 
rates in Figs. 3f and h. Even though the growth rate of C2H is 
about thirty times as efficient at vd - 30kms ' its lower abun- 
dance results in a total growth rate Tq^ lower by 0.5-1.0 orders 
of magnitude in the same region (2.5/?* < radius < 3.8/?,), 
Fig. 3i. However, at high enough drift velocities non-thermal 
sputtering is active (see below) and the net growth rate t~' is 
negative, compare with the lower-most panel in Col. 3 Fig. 2. 

The two tall peaks in Fig. 3j, between 2.5 and 3.5 /?„ are re- 
gions where the helium particles in the gas are energetic enough 
to erode dust grains. Non-thermal sputtering quickly becomes 
significant in comparison to the total growth rate Tq when the 
drift velocity reaches values around 30-35 km s"^ As such this 
process is mostly present in regions of a low dust density where 
the drift velocity may be higher. And this is probably the reason 
to why the model properties (excepting the mean drift velocity 
itself) presented in Sect. 3.3 are found to be independent of it. 
A property that is affected is the average grain radius (Fig. 3d), 
which quickly decreases by two orders of magnitude and more 
in these regions. 

That grain growth through the radical C2H is not neces- 
sarily negligible is for example seen in a check of the relative 

growth rates of model P13C16U6-J In particular there seems 

to be a tendency towards a larger importance of C2H in winds 
showing a large degree of variability in the structure. It is for 
example found to be less important in model P13C16U6-p_, 
while it is found do be the dominant growth species inwards of 
the iimermost nucleation zone in model P13C16U6-J_; this is 



illustrated by the solid line in Figs. 3f-h. C2H2 is, however, still 
the dominant growth species in the improved drift models. Note 
that our use of an equilibrium gas chemistry may result in in- 
correct abundance ratios between different molecules. Reliable 
quantitative estimates of the relative contributions of C2H and 
C2H2 to the growth rates can therefore not be made at present. 

5. Conclusions 

An important part of recent models of AGB star winds is a 
time-dependent formulation that can account for, e.g., dust for- 
mation occurring in non-equilibrium, and formation and prop- 
agation of shock waves. To the authors' knowledge there has 
not been any time-dependent wind model that self-consistently 
treats grain drift in the dust formation process. This study is 
intended to fill this gap. We have covered both physical issues 
and numerical modeling, including a discussion on the impli- 
cations and relevance for the wind structure. The new models 
are based on the wind model descriptions introduced by Sandin 
& Hofner (2003a,b, Paper 1 & Paper 11, respectively). 

The new models presented here have been found to re- 
produce many properties of the previous drift models dis- 
cussed in Paper 1 and Paper 11. New effects have, however, 
been introduced through the improved treatment. In particu- 
lar the micro-physical details of the dust formation process are 
changed. An example is acetylene that normally is the main 
growth species in the wind (e.g. Gail & Sedlmayr 1988). This 
molecule is only weakly bound to the surface of dust grains 
and is therefore sensitive to the drift velocity. Its role in the 
grain growth process may be taken over by other more strongly 
bound growth species when the drift velocity increases above 
about lOkms"', partly preventing a decrease in the net growth 
rate. Dust destruction by non-thermal sputtering seems to be 
insignificant in view of the generally low drift velocities. 

It should be noted that the drift-dependent dust formation 
in principle is grain size dependent through the drift velocity. 
The results that have been presented here are based on calcu- 
lations using one mean drift velocity. A more consistent treat- 
ment might result in slightly different properties of the dust and 
the wind. Such a study, however, requires a significantly larger 
computational effort and is beyond the scope of the present pa- 
per. 

The grain growth efficiency increases significantly with 
the drift velocity, even at values as low as a few kms"'. 
Consequently the improved models produce significantly larger 
amounts of dust than the previous models did. However, it 
should be noted in this context that some of the micro-physical 
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assumptions in the present models make detailed quantitative 
predictions of wind properties difficult. This concerns in par- 
ticular the assumption of chemical equilibrium in the gas phase 
and the adopted sticking coefficients which may have a notice- 
able effect on the nucleation and growth rates. Further stud- 
ies, beyond the scope of this paper, are needed to clarify these 
points. Still, in view of ffie results presented here, we conclude 
that ffie effects of drift in the grain growffi rates cannot be sim- 
ply ignored in wind models. 
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